Geometric Frustration and Dimensional Reduction at a Quantum Critical Point 
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We show that the spatial dimensionality of the quantum critical point associated with Bose-Einstein con- 
densation at T = is reduced when the underlying lattice comprises a set of layers coupled by a frustrating 
interaction. Our theoretical predictions for the critical temperature as a function of the chemical potential cor- 
respond very well with recent measurements in BaCuSi206 [S. E. Sebastian et al, Nature 411, 617 (2006)]. 
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The universality class of a critical point (CP) depends on a 
few properties such as the symmetry of the underlying model, 
the range of the interactions, the number of components of 
the order parameter (OP), and the space dimensionality df^. 
It is well known that for highly anisotropic systems such as 
weakly coupled layers, the universality class changes when 
the system approaches the CP. A dimensional crossover takes 
place: the effective dimensionality is reduced beyond a cer- 
tain distance from the CP, determined by the weak inter-layer 
interaction. Sufficiently close to the CP the transition is how- 
ever three-dimensional. In contrast to this common behavior, 
the dimensional reduction (DR) discussed in this paper occurs 
when the system approaches a gaussian quantum critical point 
(QCP). The source of this qualitative difference is in the nature 
of the interlayer coupling. We show that the interlayer cou- 
pling vanishes right at the QCP for a chemical potential tuned 
Bose-Einstein condensation (BEC) of interacting bosons. We 
then argue that this effect is relevant for the field tuned QCP 
of a geometrically frustrated quantum magnet. 

Although geometric frustration has previously been in- 
voked J2l as a mechanism for DR, zero-point fluctuations are 
expected to restore the inter-layer coupling yfl, as explicitly 
shown by Maltseva and Coleman 0]. In distinction to this ex- 
pectation, we show that this coupling is suppressed near the 
BEC-QCP, relevant to spin dimer systems in a magnetic field. 
In this case, the spatial dimensionality of the gaussian QCP is 
d = 2. Interactions between either thermally excited or quan- 
tum condensed bosons induce a crossover to d = 3 away from 
the QCP. Key to this result is the observation that zero-point 
phase fluctuations of the OP are suppressed near a chemical 
potential tuned BEC. The first experimental evidence of this 
phenomenon was found very recently by measuring critical 
exponents of a field induced QCP in BaCuSi20g JH. 

We start by presenting the rigorous result for the model of 
a chemical potential tuned BEC where the bosons are located 
on the sites of a body-centered tetragonal (bet) lattice. In the 
second part of the paper we discuss the relevance of this re- 
sults for S = h spins forming dimers on a bet lattice, closely 
approximating the case of BaCuSi206 an. This offers a 
quantitative explanation of the observed DR in this systemfl. 



We start from the Hamiltonian of interacting bosons 
H B = (£k - fi) a^ak + u mrii 



(1) 



where m = ajoi is the local number operator of the bosons 
and = ajV k ' Ri /\/iV the corresponding creation op- 
erator in momentum space. The tight binding dispersion for 
nearest neighbor boson hopping on the bec lattice is 
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where k 

£n (kill = 



I = (k x ,k y ) refers to the in plane momentum. 
= tii (2 + cos k x + cos k y ) is the in-plane dispersion 
while 7 (k||) = cos 4f cos For t», tj_ > and tn > t±/2, 
a BEC takes place at K|| = (it, ir). Since 7 (K11) = 0, the 
dispersion £k at the condensation momentum is independent 
of k z . In case of the ideal Bose gas (u = 0) this implies for 
T = that different layers decouple completely. Only ex- 
citations at finite T with in-plane momentum away from the 
condensation point can propagate in the z-direction. This be- 
havior changes as soon as one includes boson-boson interac- 
tions (u > 0). States in the Bose condensate scatter and create 
virtual excitations above the condensate that are allowed to 
propagate in the 3-direction. These excitations couple to con- 
densate states in other layers^. The condensed state of inter- 
acting bosons is then truly three dimensional, even at T = 0. 

The above argument for "dimensional restoration" due to 
interactions does not apply in case of chemical potential tuned 
BEC. In this case, the number of bosons at T = is strictly 
zero for fi < 0, i.e. before BEC sets in. The absence of par- 
ticles makes their interaction mute and one can approach the 
QCP arbitrarily closely without coherently coupling different 
layers. While the Bose condensed state for \i > and the en- 
tire regime for T > is three dimensional, the decoupling for 
(fi < 0, T = 0) has dramatic consequences. We show that the 
BEC transition temperature varies as 
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T c oc p: 2 / d holds instead for an isotropic Bose system in d > 2. 
Despite the fact that different layers are coupled at finite T the 
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BEC-transition temperature, Eq.([3j, depends on p just like the 
Berezinskii-Kosterlitz-Thouless (BKT) transition temperature 
of a two dimensional systemJH]. 

The renormalization group (RG) calculation used to obtain 
this result (a one-loop RG calculation in analogy to Ref.Jit]) 
shows that the finite temperature transition is a classical 3-d 
XY transition, not a BKT transition. We conclude, there- 
fore, that the T = QCP of chemical potential tuned BEC 
with three dimensional dispersion, Eq.©, is strictly two di- 
mensional. The system then crosses over to be three dimen- 
sional for p > or T > 0, where the density of bosons be- 
comes finite and boson-boson interactions drive the crossover 
to d = 3. The transition temperature of this three dimensional 
BEC is given by the two-dimensional result, Eq.©. It is im- 
portant to stress that the vanishing density for (p < 0, T = 0) 
implies that these results are not limited to weakly interacting 
bosonsjgl. 

The detailed derivation of Eq.(|3]i using the RG approach 
will be presented in a separate publication ifioll . Here we 
present a heuristic derivation of the same result based on an 
approach introduced by Popov 111 ill and further explored by 
Fisher and Hohenberg Q]: infrared divergencies are cut-off 
for momenta k < fco ~ ^ p/t\\. We analyze the interact- 
ing Bose system in the disordered phase and perform an ex- 
pansion in the interlayer hopping amplitude t±/tu . Dominant 
interactions at low density are given by ladder diagrams (see 
FigQJ), yielding a renormalized boson interaction (i.e. the 
scattering matrix for bosons in the same layer) 11211 : 
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for u — ► oo (hard core bosons). The bare interlayer cou- 
pling leads to scattering of bosons between different layers. 
The corresponding scattering matrices between neighboring 
layers, V\ (see FigQJ)), and second neighbor layers, v% (see 
FigOJ) are then given as (I = 1,2) 



vi~-\ — 
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(5) 



where the overall negative sign results from the fact that the 
lowest order contribution to v-y^ are of order Vq. This in- 
terlayer coupling is on the interaction level, and leads to 
new non-local interaction terms viriini + i ez in the low en- 
ergy Hamiltonian of the model. The origin of these cou- 
plings are T = quantum fluctuations of the interacting 
Bose system. Pairs of boson propagate as virtual excita- 
tions between layers and mediate the non-local boson-boson 
coupling 14]. It is crucial to observe that, no coherent boson 
fflkii ,n+( between layers emerges for T = 0. 



hopping *!,<,,„ 



Hb is invariant with respect to the discrete ^-symmetry: 
k x — > — k x + 2ir and k z — > k z + it. As long as this sym- 
metry is intact, no term t* ± 1 cos (k z ) in the dispersion is al- 
lowed, while coherent hopping between second neighbor lay- 
ers with t* ± 2 cos (2k z ) does not break the ^-symmetry. To 
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FIG. 1: (Color online)(a) Ladder diagrams that provide the dominant 
contribution to the intra-layer scattering in the low density regime 
fl2ll . (b) and (c) leading order diagrams that contribute to the coher- 
ent inter-layer hoppings tj_,i an d tj_ 2- 



determine these coherent interlayer hoppings t*\_ t we perform 
a mean field (MF) theory of the low energy problem with 
interlayer interactions i>;. We approximate viniJii + i Bz — > 



,t. 



vi \aiai + ie z ) ffljOi+i 



and obtain 
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The expectation values (a^ „ak||,n+iy for the single parti- 
cle overlap between neighboring layers are determined self- 
consistently. As expected, we find t_]_ 2 (T = 0) = and 
t_]_ 1 (T) = 0. The former result reflects the fact that no coher- 
ent motion is possible at T = 0, while the latter is caused by 
the ^-symmetry, forcing the hopping between nearest neigh- 
bor layers to vanish at all T. The solution of Eq.© for the 
coherent second neighbor hopping is 



t> 2 



til J tu tnk? 



(7) 







Using the above result for V2 it then follows t* ± 2 ~ 

T \ntJ/^ ' Since the density of bosons is p ~ 
Tin (T/fi) /tn, one sees that thermally excited bosons induce 
a coherent hopping between second neighbor layers, i.e. t\_ 2 
oc p/ In (i || I \i\ . While the amplitude of this coherent hopping 
is small, the finite T transition will be three dimensional. The 
ordering temperature of this 3-d XY transition is given by the 
MF condition: 



p c = V p, 



(8) 



and as usual for strongly anisotropic systems, its value is given 
by the characteristic temperature scale of the in plane order- 
ing. Since t\_ 2 <C en (feo), the d = 2 fluctuations dominate 
the magnitude of T c at very low densities resulting in Eq.Q. 
For the same reason, we also obtain d = 2 expressions for: 

p(T = Q,p) oc phi — 



p(T, p = 0) oc 



tii 1 



(9) 
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Based on these results we next address the origin of DR 
in the frustrated magnet BaCuSi2C>6 13]. We start from a 
Heisenberg Hamiltonian of S = | spins dimers on a bet lat- 
tice, closely approximating the case of BaCuSi20g 10, 0]. 
The dominant Heisenberg interaction, J^i S»i ■ S$2, is be- 
tween spins on the same dimer i. Since there are two low 
energy states in an applied magnetic field, the singlet and the 
Sfi + S? 2 = 1 triplet, we can describe the low energy sec- 
tor with hard-core bosons. The triplet state corresponds to an 
effective site i occupied by a boson while the singlet state is 
mapped into the empty site lfl3i [Till . The resulting low en- 
ergy effective Hamiltonian corresponds to a gas of interact- 
ing (infinite on-site repulsion) canonical bosons, as given in 
Eq.( [TJ). The number of bosons (number of triplets) equals 
the magnetization along the z-axis. The chemical potential 
fjb = gps{H — H c ) is determined by the applied magnetic 
field H and the critical field gpsH c = J — 2 J'. The hoppings 
tii = J' and t± = J 1 - are determined by the inter-dimer ex- 
change interactions between spins on the same bilayer, J' ~ 6 
K Gl [3 E3] and on adjacent bilayers, J 1 - < J'. The 
modulation of the BaCuSi2C>6 lattice structure along the c- 
axis leads to an alternation of two non-equivalent bilayers 
A and B, with intra-dimer interactions J a = 49.5(1)K and 
Jb = 54.8(1)K J3,[3]. This alternation reduces the magni- 
tude of the residual non-frustrating inter-layer couplings char- 
acteristic to all real systems Il7tl . while the principal treatment 
of BaCuSi2C>6 presented here remains unaffected. 

The correspondence between the quantum spin model for 
BaCuSi2C>6 with the boson model of Eq.(Q]i allows us to in- 
terpret T c of Eq.Q as the phase boundary as a function of 
M = gps{H — H c ). At this phase transition, we also expect 
that the Zi symmetry will be broken as well. It is interesting to 
analyze the dimensional crossover and the coupling between 
second neighbor layers directly in the spin language. For clas- 
sical spins Si at T = 0, the frustrated nature of J- 1 produces a 
perfect decoupling of the OP's (X Y staggered magnetization) 
on different layers. However, this decoupling is unstable with 
respect to quantum or thermal fluctuations [4]. Either of these 
fluctuations induces an effective inter-layer coupling via an 
order from disorder mechanism as illustrated in FigfJ] When 
the sum of the four spins on a given plaquette, Sp, is exactly 
equal to zero the coupling between that plaquette and the spins 
which are above (Sp) and below (Sp) cancels out. However, 
the effect of phase fluctuations is to produce a net total spin on 
the plaquette, Sp ^ 0. Since Sp is antiferromagnetically cou- 
pled with Sp and Sp (see Figfft)), an effective ferromagnetic 
(FM) interaction results between Sp and Sp, i.e. between 
second neighbor layers. 

The defining characteristic of the BEC quantum phase tran- 
sition is that it is driven by amplitude fluctuations of the OP, 
in contrast to the XF-like transition that is driven by phase 
fluctuations. This difference is vital to the effective coupling 
at the QCP: it vanishes due to its quadratic dependence on the 
amplitude of the OP. The remarkable consequence is a DR of 
the gaussian QCP from d — 3 to d = 2. 

Our previous analysis shows that the non-universal prefac- 



(a) 



(b) 



FIG. 2: (Color online)(a) The perfect antiferromagnetic (AF) order 
of the four spins in the square plaquette precludes an effective cou- 
pling between Sb and St ■ (b) A phase fluctuation of the AF OP 
induces an effective ferromagnetic coupling between Sb and St- 
For BaCuSi20s, each site represents a dimer. 



tors of Eqs.© and (O can be determined to high accuracy by 
using a theory for a strictly d = 2 system. The RG and MF 
approaches used to describe the quasi-condensate phase of a 
weakly interacting two-dimensional Bose gas 1 1 ill give the 
proper universality class and generic T vs. fj, dependence, yet 
are quantitatively not adequate for realistic densities. The lim- 
itation of these treatments arises from the insensitivity of the 
size of the critical region AT (where the fluctuation correc- 
tions associated with the BKT transition are important) to the 
smallness of the interaction: AT/T C ~ 1/lntii/uo. This lim- 
itation is, however, overcome by our use of Monte Carlo (MC) 
simulations to obtain the non-universal constants that appear 
in Eq.©, and the results of Prokof 'ev and Sustinov 1 1811 who 
computed these constants explicitly for realistic low densities 
and weak interactions, obtaining the following expression for 
the phase boundary of the quasi-condensate [18]: 



/'<■ 



IT J' Vq 



(10) 



where £ M = 13.2 ± 0.4. In Fig. [3^, we compare the exper- 
imental data for BaCuSi206 ||5|] with the result of Eq.(fT0l) 
and Monte Carlo (MC) simulations of hard core bosons on 
a square lattice (L x L with L = 32) with hopping t\\ = 
J' = 6K BG3]. The agreement is good for T < 200mK 
(p < 0.02) but, as expected, there is a significant deviation 
at higher temperatures (densities). The fact that the measured 
T c becomes significantly higher than the MC result at higher 
temperatures indicates that neglecting the effective inter-layer 
tunneling is no longer valid in BaCuSi206 for p > 0.02. 
Fig. shows a similar comparison for p(p, T ~ 0) and 
p(p = 0,T) [see Eqs.©]. Again, we compare the experi- 
mental data against the MC simulation because the MF ap- 
proximation that leads to Eqs.© is adequate to determine 
the generic p and T dependence, but cannot reproduce the 
non-universal constants. Our theory also predicts a linear de- 
pendence of the specific heat C(T, H c ) and the nuclear re- 
laxation time l/Ti(T, H c ) as a function of T at the QCP of 
BaCuSi 2 6 . 

To compute the exponent of the next order correction to 
Eq.© we note that the effective boson-boson interaction 
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ya(p) is obtained as an expansion in the small parameter p 1 / 2 
1 12|: v (p) = v (l + ap 1 / 2 + ...). While the first term 
in this expansion is determined by the ladder diagrams of 
FigOJ, higher order diagrams contribute to the second term. 
The MF relation ([8]) implies that the next order correction 
to Eq.dTOil is proportional to T 3 / 2 . The value of u\ deter- 
mines the crossover between the linear regime consistent with 
a d = 2-QCP and the T 3 / 2 regime characteristic of a d = 3 
BEC. Such a crossover was reported in BaCuSi 0]. By 
following a similar procedure, we can demonstrate in general 
that the phase boundary equation of a d-dimensional bosonic 
system that comprises d — 1 -dimensional regions coupled via 
a frustrated interaction is p c ~ AT^ 1 ^ 2 + BT d / 2 for low 
enough p. 
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FIG. 3: (Color online) (a) Phase boundary near the QCP measured 
in BaCuSi 2 6 H compared to the curves obtained from a MC 
simulation and Eq.JlOb for a d = 2 gas of hard-core bosons on 
a square lattice with tu = J' = 6K.(b) Similar comparison for 
p(T = 30mK lA t = H - H c ) and p(T = 30mK, H = H c ). 
We have neglected the density of bosons on the B-bilayers because 
J'' ~ 9PbH 3> \tj_ 1 1 as long as H is not close to JijgpB — 
H c + 3.4T. 

Our mapping of the spin problem to the boson model Hb 
is based on the assumption that only the lowest triplet and 
the singlet modes are important at low energies. The low 
density expansion for the boson problem is then well jus- 
tified, as the Xy-symmetry of the original spins Sj is di- 
rectly responsible for the charge conservation of Hb- Re- 
cently, it was shown by Rosch and Vojta lfl9ll that the inclu- 
sion of the two higher triplet modes generates a small coher- 
ent second neighbor hopping of low energy triplets between 
layers £j_ 2 — Jj_/J 5 . This interesting effect restores the 
d = 3 character of the spin problem. For realistic values of 
J = 49.5(1)K and J± < J', we find that J^/J 5 < O.lmK 
in BaCuSi2C>6. This implies that the mechanism discussed 
in our paper is still dominant for all experimentally accessi- 
ble temperatures T > 30mK. Moreover, the U(l)-symmetry 



breaking terms induced by dipolar interactions will produce 
a crossover to QCP with discrete symmetry at T ~ lOmK 
ll20ll before the mechanism of Ref. fl9ll sets in. Finally, the 
inevitable presence of finite non-frustrated couplings in real 
systems will eventually restore the three dimensional behav- 
ior below some characteristic temperature To (for BaCuSi206 
we estimate To < 30mK |[l7lo . We stress that our theoretical 
results for Hb [Eq.([T|l] are not affected by these considera- 
tions. There exists a non-trivial three dimensional interacting 
many body system with a strictly d = 2 QCR 

In summary, we demonstrate that the dimensionality of the 
BEC-QCP is d = 2 when the inter-layer coupling is frustrated. 
However, this coupling is relevant for changing the thermody- 
namic phase transition from BKT type to the 3d-XY univer- 
sality class. These results explain quantitatively, and with- 
out free parameters, the DR manifested in the experimentally 
measured quantum critical exponents of BaCuSi20g fl. 
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